Raman study of the phonon symmetries in BiFeOa single crystals 
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In Bismuth ferrite (BiFeOs), antiferromagnetic and ferroelectric order coexist at room temperature, making 
it of particular interest for studying magneto-electric coupling. The mutual control of magnetic and electric 
properties is very useful for a wide variety of applications. This has led to an enormous amount of research 
into the properties of BiFeOs. Nonetheless, one of the most fundamental aspects of this material, namely the 
symmetries of the lattice vibrations, remains controversial. We present a comprehensive Raman study of BiFeOs 
single crystals with the novel approach of monitoring the Raman spectra while rotating the polarization direction 
of the excitation laser. Our method results in unambiguous assignment of the phonon symmetries, and explains 
the origin of the controversy in the literature. Furthermore, it provides access to the Raman tensor elements 
enabling direct comparison with theoretical calculations. Hence, this allows the study of symmetry breaking 
and coupling mechanisms in a wide range of complex materials and may lead to a non-invasive, all-optical 
method to determine the orientation and magnitude of the ferroelectric polarization. 



Multiferroic BiFeOs (BFO) is one of the few materials that 
simultaneously exhibits a robust magnetic ordering and lar^e 
spontaneous ferroelectric polarization at room temperature| 1], 
making it of particular interest for studying magneto-electric 
coupling i^]^. The mutual control of magnetic and elec- 
tric properties is of great interest for applications in spin- 
tronics and magnetic storage media |5]. This has triggered 
significant interest in BFO, resulting in numerous studies in- 
cluding optical [v], and Raman spectroscopyfs'-^l^, the- 
oretical calculations! iH [12], thin film devices Li, 13] and 
electrical control of magnetic excitations |[l4l - [l8ll . Amongst 
these various techniques, the Raman spectrum of BFO is 
one of the most widely studied as it is a powerful tool 
to investigate phonons, magnons and their interaction (i.e. 
electromagnons).fl^- 19] Moreover, proper phonon mode as- 
signment is necessary to describe the phonons critical for the 
multiferroic behavior. However, even for measurements taken 
along the high symmetry directions of single crystals, con- 
troversy in the symmetry assignments of the phonon modes 
remains. The discrepancies have previously been ascribed to 
violation of Raman selection rules due to variations in strain 
fieldsH (i.e multidomain states) caused by polishing of the 
crystal surface. Once the symmetries are unambiguously as- 
signed, deviations in phonon mode behaviors could be used to 
detect the presence of symmetry breaking, multidomain states 
and phonon-magnon interactions. Furthermore, simply deter- 
mining the mode symmetry only allows for a qualitative com- 
parison with theoretical calculations. Whereas a quantitative 
comparison is enabled by measuring the Raman tensor ele- 
ments. 

To this end we have performed a comprehensive set of 
polarized micro-Raman spectroscopic studies of BFO single 
crystals with uniform ferroelectric polarization. Careful ex- 
amination and proper modeling of the rotational dependence 
of the Raman intensity enables us to unambiguously assign 



the (Ai, Ex and Ey) modes. Furthermore, we use the pre- 
sented model to show that slight misalignment of the crys- 
tal leads to ambiguity in the symmetry assignments. Indeed, 
our data reveal that comparison of spectra obtained for dif- 
ferent scattering geometries at a single polarization vector of 
the incoming light is not sufficient to have truly unambiguous 
mode assignment. Nonetheless, unambiguous assignment can 
be reached on the as grown single crystal when the Raman 
mode intensities as function of crystal rotation are measured 
(consistent with previous work on sapphire [20]). Hence, with 
the presented method polishing is omitted and the resulting 
ambiguity from misalignment can be avoided. 

The as-grown BFO single crystals used in this work have 
pseudocubic [100]pc facets with a ferroelectric single domain 
state|21] (see Supplemental Material 1 22]). The crystal struc- 
ture of BFO (rhombohedral distorted perovskite, R3c) shows 
a transition from high to low symmetry accompanied by the 
formation of spontaneous electric polarization below the tran- 
sition temperature Tc ~1100K|23J. The ferroelectricity is 
ascribed to lattice distortions (i.e. off-centering of the Bi-ions) 
and results from softening and subsequent freezing of the low- 
est frequency polar-phonon mode. The antiferromagnetic or- 
dering sets in below T at ~ 640 K with a large magnetic mo- 
ment of 4 /iB on the Fe-ions. Canting of the spins leads to 
a cycloidal spin structure with large period (62 nm)|23, 
rotating in the plane containing the electric polarization vec- 
tor P and cycloid wavevector q. At room temperature BFO 
has a perovskite pseudocubic unit cell (a ~ 3.96 A) elongated 
along the (1 1 l)pc direction coinciding with P. The point group 
is Cs^, with 13 Raman active modes, of which four have Ai 
symmetry (i.e propagate along the c-axis) and nine have either 
Ex or Ey symmetry (i.e. propagate in the x-y plane), which 
are doubly degenerate. When the laser is not along the c-axis, 
phonons can propagate in the x-z plane, which could lead to 
LO-TO splitting (i.e. lifts the degeneracy) and hence, the pres- 
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ence of A (TQ) m odes in the XX and E(LO) modes in the XY 
geometry ll25l427n . which further compHcates the analysis. 

The Raman spectra were taken in a backscattering config- 
uration with a Horiba Jobin Yvon LabRam microscope with 
a 532 nm excitation source and a lOOx objective (0.8 NA), 
resulting in a collection area of ~1 jim (see Supplemental 
Material [22]). All data presented in this work are taken at 
room temperature. Furthermore, we investigate the polariza- 
tion dependence of the Raman spectra by linearly polarizing 
the excitation laser in the plane of the sample and rotating the 
polarization direction with steps of 10 degrees over a total of 
180 degrees. The rotation is accomplished via a A/2 Fresnel 
Rhomb and is fully equivalent to an in-plane rotation of the 
sample (see Fig. [1^ and Supplemental Material 1 22]). A sec- 
ond polarizer is used to analyze the scattered light, which is 
either parallel (XX) or perpendicular (XY) to the incoming po- 
larization direction. 
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FIG. 1 : a) Experimental setup with the Fresnel Rhomb (FR) used to rotate 
the polarization of the incoming laser. The green beam is the excitation laser 
and the red the Raman scattered light, with polarizers (P), Notch filter (BS), 
objective (O) to focus down the laser and the sample (BFO). b) Typical sin- 
gle phonon spectra in XX geometry for two different [100]pc crystals (black: 
crystal I and red: crystal II) for Raman shifts between and 650 cm~ ^ . Inset: 
full range up to 1500 cm~^. The curves are vertically translated for clarity. 



Fig. [Tt shows typical Raman spectra (XX scattering geom- 
etry) taken on two different crystals (both with a [100]pc sur- 
face). The modes below 600 cm~^ are single phonon modes 
and the broad features above 600 cm~^ (see inset FigdJ)) are 
ascribed to 2 phonon excitations, which is in agreement with 
previous reports [9] . The spectrum taken on crystal I shows 
a total of 1 1 single phonon modes (see Table [B), while crystal 
II shows a total of 13 single phonon modes (i.e. all modes 
observed in crystal I and two additional modes at 53 and 77 
cm~^, which can be seen due to the use of a better filter with 
a lower cutoff frequency). Raman intensities taken on differ- 
ent locations on one crystal and on different crystals (Fig. [TJ)) 
show similar polarization dependencies (i.e the symmetry as- 
signments are consistent), which confirms the single domain 
character of the crystals. By comparing the polarization de- 
pendence of the Raman intensities of crystal I and II we show 



how a different (but homogeneous) ferroelectric polarization 
direction influences the phonon mode behaviors (discussed in 
detail below). In Fig. [2] we show the evolution of the Raman 
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FIG. 2: The evolution of Raman spectra as the Fresnel Rhomb is rotated, a) 
XX scattering geometry, b) XY scattering geometry. The spectra were taken 
on the as-grown [100]pc surface of crystal I. 

spectra as a function of in-plane crystal rotation (i.e. rota- 
tion of the polarization direction) taken on the [100]pc sur- 
face of crystal I for the XX (FigEh) and XY (Figl2b) scatter- 
ing geometries. Furthermore, we normalize the Raman spec- 
tra at 1500 cm~^ to correct for any power fluctuations of the 
laser and for polarization dependence of the reflectivity of the 
crystal. We have also confirmed that the anisotropy of the 
optical constants lEsIl does not significantly influence polar- 
ization dependence of the Raman spectra (see Supplemental 
Material [22]). To quantitatively analyze the data the spec- 
tra are fit with multiple Lorentzian oscillators of the form: 
/(c.) = /o + E.(i 



(4(u;-£; )2+f2) ) where i is the peak number, 
/o accounts for the background, Ei is the center frequency, 
is the width, and Ai is the area of peak i. The fitting is 
done with fixed mode positions, extracting mode peak inten- 
sity (li(E^)) from the ratio between area and width of the fitted 
oscillators (i.e. I^(E^)=A^/r^). In Fig. [3] we show the mode 
intensities as determined from fitting the Raman spectra as 
function of polarization angle for three representative modes. 
The polar plots indicate the presence of exactly three differ- 
ent mode symmetries. Not surprisingly we have found that all 
modes can be sorted into one of these three types. Indeed, the 
fits in Fig. [3] show that these three types match well with what 
we expect for the A, E^^ and E^y symmetries. Here we note that 
the differences between the mode behaviors can be subtle, for 
example the XY curves for the A and the E^^ modes (see Fig. 
[3]b and d) look similar. Hence, simultaneous modelling of the 
full polarization curves for both XX and XY is necessary and 
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FIG. 3: Polar plots of the mode intensities determined from the Raman spec- 
tra (left: XX and right: XY) as function of polarization rotation for three rep- 
resentative modes (a) and b): mode @ 350 cm"-*^, c) and d): mode @ 140 
cm~^ and e) and f): mode @ 471 cm~^) measured on crystal I. The solid 
lines are fits ( 350 cm~^: A, 140 cm~^:Ey and 471 cm~^: E^) of which the 
tensor elements are indicated in Table U 




12 3 4 

Misalignment angle between 
270° Poynting vector and c-axis (deg) 

FIG. 4: a) - c) Polar plots of calculated mode intensity variations as function 
of polarization angle for A , and E^ symmetry, respectively, in the case 
of the perfect alignment of the poynting vector with c-axis (black: XX, red: 
XY). The dotted lines in b and c indicate angles for which E-modes will be 
mistakingly assigned as A-modes. (d) - f) Linear plots of the calculated in- 
tensity variations in the A, E^ and E^ modes, respectively (for XY scattering 
geometry) as function of the deviation of the c-axis from the surface normal. 
The dotted lines indicate misalignment angles (0.9 and 1.5°) for which the 
symmetry of 2 out of 3 modes will be misassigned. 



in which ol are the Cs^ Raman tensors in a trigonal basis |l3j 



only then results in unambiguous assignment of the modes. 

Moreover, we find no evidence that we are probing phonons 
that propagate in the x-z plane (i.e. oblique phonons ll26l 
[27]). Indeed, such modes would exhibit LO-TO splitting 
as seen previously, with the presence of A(TO) and E(LO) 
modes leading to reduced intensities of the A(LO) and E(TO) 
modes f27]. Hence, we observe the 13 modes expected from 
group theory and the modeling shows that the phonons trans- 
form according to the zone center modes irreducible represen- 
tations. Furthermore, on a polished c-axis surface we find the 
same number of modes as for the as grown surface with the 
modes at the same frequencies (within our resolution). 

With a closer look at the model we used for the fits in Fig. [3] 
we can explain why there is controversy in the phonon mode 
assignment in the literature. The extracted mode intensities as 
function of polarization angle are modelled using the Raman 
tensors for the C^v point group (i.e. the trigonal symmetry of 
the lattice). The Raman intensity as function of polarization 
angle can be calculated using the equationl 2QlE9h . 





(1) 



with R the matrix that rotates from cubic to the trigonal orien- 
tation, and with and e^ the polarization vectors that describe 
the scattered and incoming light, respectively. The polar plots 
for the calculated Raman mode intensity as a function of po- 
larization rotation are shown in Figs. H^-c. Here we have as- 
sumed the poynting vector is perfectly parrallel with the c-axis 
[111] of the crystal resulting in easy to distinguish behaviors 
of the modes. Indeed, in this case determination of A modes 
should be easy (Fig. |4^); their intensity is independent of the 
polarization angle and should have no measurable intensity in 
the XY geometry. However, at certain angles of crystal orien- 
tation (indicated by dotted lines in Figs. [It and c one can still 
mistake an E mode for an A mode (i.e. the mode has inten- 
sity in the XX geometry but disappears in the XY geometry). 
Additional error can come from slight misalignment of the c- 
axis with respect to the propagation of the Raman laser light 
(i.e. surface normal). Figs. |4]l-f demonstrate that, for the A, 
and E^ modes respectively, introduction of a few degrees 
misalignment can lead to large variations in the intensity of 
the phonon modes (in the XY geometry), and hence crossover 
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in mode symmetry assignments. The two dotted lines in Figs. 
|4]i-f are examples of misalignment angles (0.9 and 1.5° re- 
spectively) for which the symmetry of 2 out of 3 modes will 
be misassigned. Also, the Raman beam is typically focused 
down, resulting in an average of incident angles, which al- 
ready introduces some misalignment. Here we note that our 
method results in the same mode symmetry assignments com- 
pared to Palai et al.[9], measured on a polished c-axis surface. 
However, we also show that the standard method of just mon- 
itoring the disappearance of modes when switching from XX 
to XY scattering geometry does not provide adequate informa- 
tion to unambiguously assign the phonon mode symmetries. 
Our results not only unambiguously determine the mode sym- 
metries but also explains the controversy in the literature and 
leads to direct determination of the Raman tensor elements. 

In TableUwe show the observed phonon mode frequencies, 
Raman tensor elements (a, b, d, e and/) as obtained from 
the fits and and their symmetry assignments for the [100]pc 
surface of crystal I. We provide the corresponding polariza- 
tion curves in the Supplemental Material [22]. Moreover, on 
crystal II, 13 modes were observed (i.e. the correct amount 
according to group theory) at the same frequencies and with 
the same assignments as presented in Table IH two additional 
modes were observed at 53 and 77 cm~^, the tensor elements 
and symmetry assignments for all modes observed on crys- 
tal II are shown in Table TI in the Supplemental Material [22"]. 
Here we note that the symmetry assignment of the mode at 
279 cm~^ remains challenging, because it is very weak and 
shouldering the very strong Ey mode at 288 cm~^.We have 
also checked the Raman spectra on a polished [1 1 l]pc (i-e. c- 
axis) surface (data not shown) on which we observed a total 
of 13 modes (the mode at 53 cm~^ disappeared while an ad- 
ditional mode appeared at 70 cm~^). The mode at 70 cm~^ 
remains unassigned, it probably also exists on the [100]pc sur- 
face but is too weak and close to a strong E-mode to be clearly 
visible. Furthermore, it is possible that the mode at 53 cm~^ is 
indicative of a violation of Raman selection rules due to sym- 
metry breaking. Modes at this Raman shift have been previ- 
ously assigned as A(TO) modes [9], however they should not 
be visible in our scattering geometry and we do not see evi- 
dence of the other A(TO) modes in our spectra. Alternatively 
this mode may be an electromagnon[17J. Future low temper- 
ature studies, where the linewidths are narrow, would help to 
better assign these modes. Nonetheless, using the presented 
method we have unambigously assigned the phonon modes 
and extracted the Raman tensor elements providing quantita- 
tive information for direct comparison with theoretical pred- 
ications. Furthermore, the ratio between the Raman tensor 
elements a and b are identical for the A-modes observed on 
both crystals. However, we do observe some differences be- 
tween the ratios of the E-mode tensor elements between the 
measurements taken on crystal I and II, which in no way in- 
fluences the consistency of the symmetry assignments. These 
differences may indicate that the two crystals (both are sin- 
gle domain) have a different direction and/or magnitude of 
the ferroelectric polarization. This would indeed affect the 



E-modes but not the A-modes, since the A-modes are fully 
symmetrical and constitute vibrations along the c-axis (i.e. 
parallel to the ferroelectric polarization direction). However, 
this could mean that changes in the ferroelectric polarization 
direction leave the mode symmetries unaltered. Hence, one 
needs the method presented here to observe this subtle effect 
(i.e. changes in the tensor element ratios of the E-modes) of 
different ferroelectric polarization on the Raman intensities. 
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TABLE I: Phonon mode frequencies, the Raman tensor elements (a, b, d, e 
and/) for the modes as obtained from the fits and the symmetry assignments 
for crystal I. The data for all the modes are presented in the Supplemental 
Material |22]. 

We have measured the evolution of polarized Raman spec- 
tra of BFO single crystals and extracted the polarization 
curves for every single phonon mode for both the XX and 
the XY scattering geometry. We fit the XX and XY curves 
simultaneously for each mode using a model based on the 
Raman tensors of the Csy point group (eq. [T]). As a result 
unambiguous symmetry assignment and determination of Ra- 
man tensor elements of the phonon modes is accomplished 
even on the as grown [1 0]pc surface. In Fig. [5] the ex- 
cellent and unambiguous agreement between the experimen- 
tal results and our calculations for a [100]pc surface demon- 
strate the importance of performing the Raman measurements 
over a full rotation of the crystal. Whereas the calculations 
in Fig. m demonstrate that only measuring the XX and XY 
spectra on a [1 1 l]pc surface for a single polarization direction 
(as is typically done), can easily lead to wrong assignments 
of the phonons due to misalignment of the crystal. It is clear 
that unambiguous mode assignment can only be reached if 
one monitors the Raman signal as function of rotation of the 
crystal. Simply comparing XX and XY scattering geometries 
for one polarization angle is not enough even for a c-axis sur- 
face. Besides obtaining unambiguous mode assignment for 
BFO, this work has wider implications as well. The method 
can be used on any material to check crystal symmetry and 
assign the phonon modes. Furthermore, once unambiguous 
assignment has been accomplished one can use the presented 
method to investigate symmetry breaking as well, for example 
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by studying deviations in the tensor element ratios, the sym- 
metry assignments and through observation of more than the 
predicted number of modes. This gives us a powerful tool 
to investigate occurrence of (electro)magnons and compare to 
existing reports[17, 18]. Moreover it would allow study of 
coupling mechanisms in complex materials such as multifer- 
roics, as well as provides quantitative information for direct 
comparison with theoretical predications. 
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